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Abstract 

We have analyzed the two-dimensional antiferromagnetic Heisenberg model 
on the triangular lattice using a Schwinger boson mean-field theory. By ex- 
panding around a state with local 120° order, we obtain, in the limit of infinite 
spin, results for the excitation spectrum in complete agreement with linear 
spin wave theory (LSWT). In contrast to LSWT, however, the modes at the 
ordering wave vectors acquire a mass for finite spin. We discuss the origin of 
this effect. 
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INTRODUCTION 



The two-dimensional spin-| Heisenberg antiferromagnet (HAFM) has received much at- 
tention during the last few years. The main reason for this new interest stems from the 
discovery of the ceramic superconductors, which has led to extensive studies of the HAFM 
on a square lattice [|l|. An important focal point has been the experimentally observed com- 
petition between superconductivity and antiferromagnetic order in these materials. This has 
encouraged numerous investigations of disordering effects, in particular the role of frustrating 
interactions. 

To elucidate the role of frustration, it is interesting to study cases where frustration is 
intrinsic to the lattice structure such as the HAFM on a triangular or Kagome lattice. Early 
papers on the triangular HAFM 0|| suggested that the spin-^ ground state is a spin liquid, 
lacking long-range order (LRO). More recent studies HHHJ^'i'll point to the scenario of 
an ordered ground state possibly quite close to a disordering transition. Yet other studies 



support the original proposal of a disordered ground state [10,11 



The spin dynamics of the model, intimately linked to the character of the ground state, 
also remains an open question. The problem of the spin dynamics above a ground state 
with Neel-type order (a 120° state) has been attacked with a variety of methods, like linear 
spin- wave theory (LSWT) []5J, nonlinear a model |T2j and Schwinger-boson mean-field theory 
(SBMFT) [|Bi[l3|. The SBMFT has proved successful in incorporating quantum fluctuations 
|14| , but the choice of mean-field parameters is a delicate matter |T5|JT6| . Yoshioka and 
Miyazaki || have used the SBMFT combined with a resonant valence bond ansatz, and a 
similar theory including a ferromagnetic component in the mean-field expansion is recently 
analyzed by Lefmann and Hedegard |L3]]. A drawback of these applications of SBMFT 
is that the LSWT structure of the dispersion relation cannot be recovered in the limit of 
infinite spin. Specifically, only two massless spin wave modes can be obtained, in contrast 
to the three massless modes of LSWT. This is a weakness, since, at an analytical level, the 
only reliable way to motivate the validity of the theory is to use an approach that correctly 
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reproduces the limit of classical magnetism [16 . 



We here employ another mean-field treatment, where we use a spin-rotation to identify 
the relevant mean-field parameters, treating the spin dynamics above a state close to the 
classical ground state. In this way three modes are obtained, two of which exhibit a massgap 
at finite spin. However, in the limit of infinite spin, the gap asymptotically vanishes on the 
relevant energy scale and the corresponding spin wave velocities exactly approach those of 
LSWT. We discuss the origin of massgaps within this version of SBMFT and "explain" why 
there are no gaps in LSWT at any value of the spin. 

THE MODEL 

We consider the Hamiltonian 

T~t = Jl ^2 ■ Sfi+a , (1) 
Ret 

defined on a triangular lattice and with J\ > 0. The sums run over all lattice sites R and 
three of the six nearest-neighbor vectors a. We choose the vectors a to be 



1 _ \/3 
2 Gx ~ ~2 
1 \/3 



oli = ~-e x —e y , (2a) 



a 2 = --e x + —e y , (2b) 
«3 = e x , (2c) 

in units where the lattice constant equals 1 (see Fig. p. 

The triangular lattice is a Bravais lattice and can be divided into three equal interlacing 
sublattices (tripartite lattice). In Fig. [l| the three sublattices are indicated. It follows that 
in the classical limit of large spin, the ground state is a 120° state, i.e. a coplanar state with 
spins on each sublattice being ferromagnetically ordered and rotated 120° vis a vis the spins 
on the two other sublattices. (The lattice constant of a sublattice is a/3 times the nearest 
neighbor spacing why this state also is called a "a/3 x a/3 state".) 
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When expanding around such a 120° state, it is convenient to start with a local rotation 
of the spin coordinate system. The spins on neighboring sublattices are rotated by +120° 
around the y-axis when going from R to R + a. 0: 

S R ^cos{Q-R)S R + sm{Q-R)S z R , (3a) 
S R - S R , (3b) 
S Z R -»• cos(Q • R)S Z R - sin(Q • R)S R , (3c) 

where Q is defined by Q ■ ex. = ^ + 27m with n integer. This implies that 

So v / qz Q2 | f qz qz I ni \ , rtj/ qJ/ 

R - 'Jfi+a — ► n [•JR'JR+a. ~ ^ R^ R+ct ) ~ n [^R^R+a D R D R+cJ ' ^R^R+a V 1 ) 



2 \^R^-R+a u R u R+a) ^ 

Let us pick a specific 120° state, with the (unrotated) spins on the first sublattice pointing 
in the z-direction while on the other two sublattices they span the x — z plane with directions 
— 120° (+120°) with respect to the z-axis when going in the ct (—ct) direction in the lattice. 
Any classical ground state can clearly be reached from this state by a (global) 5*0(3) rotation 
of the spins. We immediately see that the spin rotation in Eqs. (^) transfers this 120° state 
into a ferromagnetic state with quantization axes lying in the x — z plane. The first term in 
Eq. (Hj) vanishes when acting on a ferromagnetic state, and we may hence neglect it when 
modeling states with a local structure close to a 120° state. In the rotated basis, the effective 
Hamiltonian is thus given by 

^ = ^i E ~o (Sr.Sr.+<x + SrSr+o) + SrSr+c* ■ (5) 



We have verified that in our case this approach is equivalent to that in Ref. [17] where 
an analysis of the frustrated model on the square lattice is performed. It is interesting to 
note that in LSWT the term that corresponds to ^ (j^R^R+a ~ SrSr+o) m Eq- (f§) is also 
neglected since it is cubic in the Holstein-Primakoff bosons. 



To go in the +a and in the — a directions is not equivalent since the 120° state breaks parity. 
Note that in two dimensions a parity transformation is defined by a reflection in a line. 
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The Schwinger-boson calculation that follows is similar to the treatment of the J\ — J 2 
model on the square lattice |jTB[ and on the honeycomb lattice |TjJ. The spin operators Sr 
at each lattice site are replaced by two species of Schwinger bosons V R (/x = 1, 2) via 

S R = 2 6 U a mAk > ( 6 ) 

with the local constraints oi R b R = 2S. Here er = (a x , a y ,a z ) is the vector of Pauli matrices, 
and summation over repeated (Greek) indices is implied. This leads us to the following 
Hamiltonian: 

n = -JiJ2 (-|><« + iw|, a ) + £ Ar[&U 6 ^ - 25] , (7) 

where we have included the local constraints with Lagrange multipliers Xr at each site. The 
expression for the summands is 

^R, a = \ '■ X Ra X Ra '■ -S 2 , (8) 
with X Ra any of the two link operators 

^R a = b iR b iR+ a + b 2 Rb 2 R +a , (9a) 

Br<x = ArWr+c + &2H&2JI+« • ( 9b ) 

The mean-field theory is finally generated by the Hartree-Fock decoupling 

: X R,a X Ra : ~* X Ra( X Rcx) + (^Ro^^Ra ~ ( X Ra)( X Ra) > 

where the link fields (mean-field parameters) Q\ = (A Ra ) and Q2 = (B Ra ) are taken to be 
uniform and real. In our mean-field treatment, we also replace the local Lagrange multipliers 
^r by a single parameter A. 

For comparison of our results with those of LSWT || it is convenient to Fourier-transform 
the Schwinger bosons independently on each sublattice, 

6 mR» = "7^ S e~ ik Rma m» k , m = 1, 2, 3 , (11) 
y-N/2 k 



with Ri, R 2 and R 3 on the three different sublattices respectively. We then need to decouple 
the different bosons on the sublattices. This is done by the canonical transformations 

a ink = ^= (Pi,* + + P^k) , ( 12a ) 
a 2 M fc = ^ (Pi»k + 3* Pi»k + 3 P^k) , (12b) 
a 3 M fc = ^ (A^fc + 3 P2^k + 3* P^k) , (12c) 



where j = exp(i- 3 

Finally we use the Bogoliubov transformations 



Piixk = cosh ^ifcCi Mfc + sinh^i fe cl M _ fe , (13a) 
^fc = cosh ^2fcC 2Mfc + sinh^ 2fc 4 M _ fc , (13b) 
P^k = cosh^ 3fc c 3Mfc + sinh^ 3fe c 2M _ fe , (13c) 



with 



tanh(2tf mfc ) = - l Jl 3 Q T ^ m ' k , m = 1,2,3 , (14) 

to diagonalize the resulting Hamiltonian. This yields free bosons c m ^ k with dispersion rela- 
tions 

u mk = \/(A + | JiQ2lm,k) 2 ~ (|^iQi7m,fc) 2 , m=l,2,3 , (15) 
where the geometrical factors are given by 

71, fc = lj2 cos ( k ' a ) ' ( 16a ) 

1 27T 

72, fc= -^cos(fc-a + — ) , (16b) 
7s,fc= ^^cos(fc-a- y) . (16c) 

By minimizing the free energy we obtain the three equations needed to determine the 
mean-field parameters: 
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1 r d 2 k 

g XI y ^4" cosh(2$ mfc )(n mfc + |) - (S 1 + §) 







(17a) 



:E ~T sinh ( 2 ^mfc)7m,fc(n mfc + \) ~ \Ql 



m=l 








(17b) 




2 / -j- cosh(2$ mfc )7 mifc (n mfc + ~) - |Q 2 







(17c) 



Here A = 8tt 2 /(3\/3) is the 8ie£i of the reciprocal unit cell of a sublattice, and fi m |- 



[exp(/3u; m fe) — 1] _1 is the Bose occupation number. Solving the mean-field equations numeri- 
cally yields values for Qi, Q 2 , and A, which are used to determine thermodynamic quantities 
at finite temperatures. 

Given a Bravais lattice, as in the present division into sublattices is not neces- 

sary when Fourier-transforming (as it is on a non-Bravais lattice) and the theory may be 
formulated in terms of only one set of bosons. We have chosen to transform each sublattice 
independently in order to make the comparison with LSWT results more transparent. In 
our treatment we thus have a small Brillouin zone and three sets of bosons (picture A), 
which is equivalent to picture B where the larger Brillouin zone of the real lattice and one of 
the sets of bosons (say c 1 k ) is used (see Fig. ^]). The equivalence of the two pictures makes 
it evident why c 2 ^ k has to be combined with c 3 ^_ k in Eqs. ( |13b| ) and ( |13c| ). It also explains 
the structure of the geometrical factors in Eq. ( |TED (remembering that Q • a = ^). Until 
now we have used picture A, but e. g. in the figures we will use picture B (one mode) instead 
of showing plots for three different modes. Since mode 2 and 3 in picture A are degenerate 
we will sometimes refer to them as the corner modes which is picture B language. To avoid 
confusion we indicate in the following which picture we use. 



In two dimensions there can be LRO only at zero temperature. This means that exactly 
at T = there is an abrupt phase transition, and we must be careful when reducing T to 
zero. As we see from Eq. (ji7|), the parameters Qi, Q2 and A will depend on temperature 



ZERO TEMPERATURE FORMALISM 
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via the Bose occupation number n m fc, and hence, to calculate properties at T = we need 
to pass to a zero temperature formalism. Following Ref. [[2(J we note that LRO corresponds 
to a condensation of the Schwinger bosons at some points K m in the Brillouin zone and we 
obtain the mean-field equations at T = from Eq. flTT| ) by the replacement 

^^s* m 5^(k-K m ) , (18) 

^>mk 

where is a new unknown quantity measuring the Bose condensate. Since the condensation 
of the bosons implies that there is no gap in the spin- wave spectrum {ui m k=K m = in fll5l) ) at 
T = 0, these m equations together with the three equations in (|17]) are sufficient to determine 
the 3 + m parameters A, Q%, Q 2 and S^. The procedure is the following: We start at a finite 
temperature and use Eqs. (|17D to determine the mean-field parameters which gives us the 
dispersion relation. We then lower the temperature and examine the dispersion relation to 
identify the points in the Brillouin zone where the gap scales down with temperature. In 
our case this happens for each of the three modes (picture A) at k — 0. We expand the 
dispersion relation around these points and for small k the dispersion relations (|15|) takes 
the relativistic form u m k = c m \j (M m c m ) 2 + \k\ 2 . Carrying out the expansion we find that 
the spin-wave velocities c m are given by 

ci = v / (!^i) 2 |-|^ 2 (A + |J 1 Q 2 ) , (19) 

and 

c 2 = c 3 = v / (I^i)H + ^ J i^(A-|Jig 2 ) . (20) 

The masses M m in the energy gap A m = M m c^ of the spin-wave excitations can in turn be 
extracted from 

A: = v^A + f</iQ 2 ) 2 - (f^iQi) 2 , (21) 

and 

A 2 = A 3 = v / (A-|Jig 2 ) 2 -(|JiQi) 2 . (22) 
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To examine whether we have condensation of bosons for all modes, we first consider the 
gaps, or the masses, since condensation of the bosons at K m = implies that we have no 
gap in the spectrum at this point. To get all modes massless (Ai = A 2 = A 3 = 0) we need 
to have Q% = Q 2 = § (from Eqs. ( [21] ) and (|22|)). Since 7i >0 = \ and 72,0 = 73,0 = — \ we 
see from Eq. (|T4| ) that we must have equalities between finite quantities as u;ifcsinh(2'#io) = 
u lk cosh(2?9io) and a; 2 fcSinh(2i9 2 o) = ^3fc sinh(2^ 30 ) = -io 2k cosh(2$ 20 ) = -cj 3fc cosh(2^ 30 ). 
But doing the replacement in Eq. (|18D into Eq. (|i~7f) for all three modes then gives Qi ^ Q 2 
from Eqs. (|17b|) and ( |17c|) , in contradiction to the initial assumption. The only way to get 
three massless modes is to have uik sinh(2-(9io) = ujik cosh(2^io) infinite, which then yields 
Qi = Q2 at infinite spin S as seen from Eq. ( |17a|) . 

The result of the procedure above is that for finite spin S we will get a condensation 
of the bosons only for the first mode which scales faster to zero than the two other modes. 
For this mode we do the replacement shown in Eq. ([18]) (m = 1 and K\ = 0) into Eq. (|17D 
(5*2 = S£ = 0) and since A x = 0, Eq. (^) gives A in terms of Qi and Q 2 . After these 
substitutions we have no temperature dependence in Eqs. (|iT|). 

With this zero temperature formalism we have studied the spin dependence of the mass 
M 2 (= M 3 ). As we see in Fig. [5] the mass goes to zero in the limit of infinitely large spin 
and we recover the three massless modes of LSWT. The dispersion relation in Eq. ([15|) can 
be rewritten as 

^ = vU + §(3Qi + Q2hm, k )(j- 1 - §(3Qi - Q 2 h m , k ) , (23) 

which immediately can be compared with the LSWT result, Eq. (13) in Ref. 0. To get the 
LSWT results we need Q\ = Q 2 = which is (as we have already seen) the condition 
to obtain three massless modes (Ai = A 2 = A 3 = 0). Thus, in the limit of infinitely large 
spin, where we in fact obtain three massless modes, we exactly recover the LSWT results. 
In Fig. |] we show the dispersion relation for two values of the spin S. We see that on the 
relevant energy scale the gaps at the corners (picture B) go to zero and for infinite spin we 
recover the LSWT dispersion relation. 
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We have also calculated the spin-wave velocities c m (m = 1, 2, 3) at T = 0. In Fig. |5] we 
see the results compared with the LSWT results extracted from Ref. || (which agree with the 
nonlinear a model results in Ref. |TJ): c\ SWT = ^J x 5 and 4 SWT = C p WT = c\ SWT /V2. 
For large spin S we have good agreement with these results. We have also examined the 
ratio of the two spin-wave velocities, as shown in Fig. |6|. 



DISCUSSION 

To summarize: We have performed a SBMFT analysis of the spin dynamics of the 
antiferromagnetic Heisenberg model on a triangular lattice, expanding around a state with 
local 120° order. Considering that the 120° state has been favored as a candidate ground 
state in several studies, this is a natural approach. In order to identify the relevant fields in 
such an expansion, we have performed a spin-rotation and we have also neglected the part 
of the Hamiltonian which gives no contribution in the classical limit. 

In the limit of infinitely large spin we obtain three massless modes in the dispersion 
relation of the spin-waves, exactly recovering the LSWT results, including the spin wave 
velocities. This is an improvement on earlier SBMFT calculations on the triangular lat- 
tice |||13], where the S — > oo limit does not come out correctly. At finite spin, however, 



two of the modes acquire a mass. This can be understood as follows: The effective Hamilto- 
nian in Eq. ([5]) is invariant under rotations in the x — z plane, and the SU(2) symmetry of 
the original model is lowered to U(l). Since the 120° state breaks this U(l) symmetry, we 
should obtain one Goldstone mode, which agrees with our results for finite spin. We do not 
believe that the massiveness of the other modes are due to enhanced quantum fluctuations 
at finite spin: If the effective Hamiltonian had preserved the SU(2) symmetry of the full 
Hamiltonian we should probably have recovered the three Goldstone modes as expected. We 
are also led to this conclusion by the discussion in Ref. where the analogous problem of 
the frustrated square lattice model is considered. 

We have already mentioned that the term in the full Hamiltonian that we neglect in 
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order to obtain the effective Hamiltonian is also neglected in LSWT. Yet, LSWT predicts 
three massless modes for arbitrary spin, even though the LSWT Hamiltonian is not SU(2) 
invariant. From Fig. |3] we see that for large spin the mass (= M3) goes approximately 
as S~ 1 ' 5 and from Fig. |5| we see that the spin-wave velocities goes like the LSWT result, 
that is, they are proportional to S. Thus, the energy gap goes like S"" 1 ' 5 * S 2 = S* 5 . Since 
the energy in LSWT only retains powers of 5" larger than or equal to 1 and the massgap is 
approximately proportional to S ' 5 , one still obtains massless modes in LSWT. In the same 
sense that LSWT is correct for large spin (neglecting terms of order S n , n < 1) our results 
are valid for large spins, neglecting energy terms of the order of the massgap. 

It is interesting to note that the modes at the corners of the Brillouin zone (picture B) 
describes fluctuations in the ^-component of the spin, while the center mode describes in- 



plane fluctuations fl2| . The first term in Eq. (y) is independent of the ^/-component, which 
means that fluctuations in this component are described equally well by the effective and 
the full Hamiltonian. This makes the value of the spin-wave velocity at the corner modes 
more reliable than that at the center mode, since this last value directly depends on the 
difference between the effective and the full Hamiltonian. 

It should be interesting to perform an extended analysis of the spin dynamics above 



a state with local 120° order. Starting with the approach of Lefmann and Hedegard |L3 



one may write the Hamiltonian in terms of a singlet field (also considered by Yoshioka and 
Miyazaki ||) and a field representing the ferromagnetic correlations, and then perform a 
spin rotation. Certain pieces of the two fields then combines to a constant term and the 
remaining parts can be written in terms of four effective fields. Diagonalizing this form 
of the Hamiltonian using the SBMFT may possibly yield a dispersion relation which can 
take all solutions into account within the same mean-field theory. The solution presented in 
this paper is valid for large spin, as we have already argued, but this approach may help 
answering for which magnitudes of the spin the solutions of Ref. || and Ref. |13[ are relevant. 
We will return to this subject. 
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FIGURES 

FIG. 1. Triangular lattice. The vectors a point to half the number of nearest neighbors of a 
site. Sites of the three sublattices are indicated by squares, crosses and circles respectively. 

FIG. 2. The Brillouin zones of the triangular lattice (large) and of the sublattices (small). The 
numbers 1 to 3 are indicating the parts of the large Brillouin zone that corresponds to the three 
different modes in the small Brillouin zone. Three possible vectors Q (see text) are also shown. 

FIG. 3. The mass M2 (= M3) of the corner modes (picture B (see text)) vs. spin S. 

FIG. 4. The dispersion relation for spin a) S = 0.5 and b) S = 200 at zero temperature. The 
gap at the "ferromagnetic" point is always zero while the gaps at the " antiferromagnetic" points 
only tends to zero for large spin. We use picture B (see text). 

FIG. 5. The ratio of the spin-wave velocities c m m = 1, 2, 3 obtained by our theory (SBMFT) 
and by LSWT (or by Ref. @) vs. spin S. 

FIG. 6. The ratio between the center mode spin- wave velocity c\ and the corner mode (picture 
B (see text)) spin-wave velocities C2 (= C3) vs. spin S compared with the LSWT result y/2. 
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